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Abstract. We construct an explicit solution of the Cauchy initial value problem for the n-dimen- 
sional Schrodinger equation with certain time-dependent Hamiltonian operator of a modified oscil- 
lator. The dynamical ST/ (1,1) symmetry of the harmonic oscillator wave functions, Bargmann's 
functions for the discrete positive series of the irreducible representations of this group, the Fourier 
integral of a weighted product of the Meixner-Pollaczek polynomials, a Hankel-type integral trans- 
form and the hyperspherical harmonics are utilized in order to derive the corresponding Green 
function. It is then generalized to a case of the forced modified oscillator. The propagators for 
two models of the relativistic oscillator are also found. An expansion formula of a plane wave in 
terms of the hyperspherical harmonics and solution of certain infinite system of ordinary differential 
equations are derived as a by-product. 



1. Introduction 

The time-dependent Schrodinger equation for a free particle 

iipt + Atlj = 

in the Euchdean space of n dimensions R"" can be rewritten in a Hamiltonian form as 



n 



s=l 



where aj and are the creation and annihilation operators, respectively, given by 
as in [26]. They satisfy the familiar commutation relations 



al, al, 



0, 



[as, as'] 

which are invariant under the transformation 

as as (t) 



as, al^ 



Sss' 



The substitution 



(s,s' = 1,2,... ,n) 
I - at (t) = ate-*. 
H^H{t) = -J2{<^sit) + al{t)Y 



e'^as 



s=l 



(1.1) 

(1.2) 

(1.3) 

(1.4) 
(1.5) 
(1.6) 
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results in the time-dependent Schrodinger equation for a modified oscillator 

i^-± = H{t)i^ (1.7) 

with the Hamiltonian of the form 

s=l s=l s=l 

In the present paper we construct an exact solution of this equation subject to the initial condition 

i^{x,t)\,^^ = ijo{x), (1.9) 

where ipo (x) in an arbitrary square integrable complex valued function from {R"") ■ The explicit 
form of equation (11.71) is given by (16.11) below. 

The paper is organized as follows. In section 2 we remind the reader about the solution of 
the stationary Schrodinger equation for the n-dimentional harmonic oscillator in hyperspherical 
coordinates and discuss the corresponding dynamical SU (1, 1) symmetry group. In section 3 we 
consider the eigenfunction expansion for the time-dependent Schrodinger equation (11.71) . The series 
solution of the initial value problem (ll.7p - (ll.9l) is obtained in section 6 after discussion of the 
Meixner-Pollaczek polynomials and evaluation of the Fourier integral of their weighted product in 
sections 4 and 5 respectively. Finally we construct the corresponding Green function in section 8, 
while introducing a required integral transform in section 7. A generalization to the case of the 
forced modified oscillator is given in section 9. An expansion formula for the plane wave in terms of 
the hyperspherical harmonics in and its special cases are discussed in section 10. A certain type 
of the time-dependent Schrodinger equation is considered in section 11 in an abstract form. The 
propagators for two models of the relativistic oscillator are derived in the next two sections. An 
axillary solution of an infinite system of the ordinary differential equations is found in section 14 as 
a by-product. Appendix at the end of the paper contains another required integral evaluation. 

The exact solution of the one- dimensional time- dependent Schrodinger equation for a forced 
harmonic oscillator is constructed in [22], [23], [21], and [3S]; see also references therein. These 
simple exactly solvable models may be of interest in a general treatment of the non-linear time- 
dependent Schrodinger equation; see [32], [3S], [IS]' [52], [51], [13] and references therein. They may 
also be useful as test solutions for numerical methods of solving the time-dependent Schrodinger 
equation. 

2. Dynamical Symmetry of the Harmonic Oscillator in n- Dimensions 

Our time-dependent Hamiltonian operator (11.81) has the following structure 

H{t) = Ho + H,{t), (2.1) 

where 

1 " 

Ho = -J2{asal + ala,) (2.2) 

s=l 

is the Hamiltonian of the n-dimensional harmonic oscillator and 

s=l s=l 
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is the part depending on time t. 

The stationary Schrodinger equation for the harmonic oscillator in n-dimensions 

92 



Ho^ = E^, Ho = IJ2(-7^ + ^"] (2-4) 



can be solved explicitly in Cartesian and (hyper) spherical coordinate systems; see, for example, 
[6l. [55] and references therein. 



In spherical coordinates r,Q given by a certain binary tree T, see [16], [55], [61] and references 
therein for a graphical approach of Vilenkin, Kuznetsov and Smorodinskii to the theory of (hy- 
per)spherical harmonics, we look for solution in the form 

ij = YKA^) R{r), (2.5) 

where are the spherical harmonics constructed by the given tree T, the integer number K 
corresponds to the constant of separation of the variables at the root of T (denoted by K due to 
the tradition of the method of i^-harmonics in nuclear physics [S3]) and u = {hyh, ■■■ , ^p} is the 
set of all other subscripts corresponding to the remaining vertexes of the binary tree T. The radial 
wave function R (r) , which satisfies the normalization condition 

POO 

/ R^ (r) r"-Mr = 1, (2.6) 
Jo 

is as follows 



where {^) are the Laguerre polynomials; see [T], [2], [B], [H], [21], 
[6^ . and references therein for the advanced theory of the classical orthogonal polynomials. 



The corresponding energy levels are 

E = N + n/2, (N - K) /2 = k = 0,l,2,... (2.8) 

and the normalized wave functions are given by 

* = "i^NKu (r, n) = Yku (n) Rnk (r) , (2.9) 

where Yku (^) are the spherical harmonics associated with the tree T and the radial functions 
Rnk {f^) are defined by (12.71) . The wave functions of the one- dimensional harmonic oscillator 

(x) = , ^ e-^"'''HN (x) (2.10) 

can be obtain from 02.91) by letting n = 1 and i^' = 0, 1 and invoking the familiar relations 

H,, (0 = (-1)' 2''k\ L,"' {e) , H,u^, (0 = (-1)^= 2''^'k\ CLf if) (2.11) 
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between the Laguerre (^) and Hermite (0 polynomials, respectively. Thus 

Rnq {x) , 



^2 -^^1 ' 

for even and odd A^, respectively; see p6] and for more details. 



(2.12) 



The n-dimensional oscillator wave functions (12.91) have the following group-theoretical properties. 
Introducing operators 



^+ = ^E J- = \Y.^^s)\ (2.13) 

s=l s=l 
^ n ^ 

■^^o = 2 5Z (^^^^ + ^^^^) = 2^°' 

s=l 

one can easily verify the following commutation relations 

[Jo, J±] = ±J±, [J+, J-] = -2Jo. (2.14) 

For the Hermitian operators 

J. = l{J+ + J-), Jy = j^{J+~J.), J-. = Jo, (2.15) 

we get 

\Jxi Jy\ '^Jzi \'^yi z\ '^Jxi \Jzi Jx\ iJy (2.16) 

These commutation rules are valid for the infinitesimal operators of the non-compact group SU (1,1); 
see, for example, [13], [25], [16] and [55] for more details. 

We are going to use a different notation for the wave function (12.91) as follows 

V'jm {x) = i'j-m{u} (x) = ^NK,y (r, n) = Yku {^) RnK {t) , (2.17) 

where the new quantum numbers are j = K/2 + n/4 — 1 and m = N/2 + n/A with m = j + 1, 
j + 2, ... . The inequality m > j + 1 holds because of the quantization rule (12. Sp . which gives N = K, 
K + 2,K + 4,.... 

The operators J± and Jo in the spherical coordinates r, Q have the form 

1 / n d \ 1 

J±=2(ffo-r^±2±rg^), J„ = -^H„ (2.18) 

and their actions on the oscillator wave functions are 

J±^im = V (m =F j) {m±j± 1) 1pj,m±l, Joi^jm = mipjm, (2.19) 

whence 

^Vjm = j(j + l)V'im (2.20) 

with J^ = Jg + Jo — J-J+ = Jq — Jo — J-^-J- . These relations coincide with the formulas that 
define the action of the infinitesimal operators J± and Jo of the group SU {1, 1) on a basis \j,m) 
of the irreducible representation D;^ belonging to the discrete positive series in an abstract Hilbert 
space [13]. In our realization of this basis in terms of the wave functions (12.171) . depending on 
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the number n = dimii", the moment j = K/2 + n/A — 1 of the group SU (1, 1) and its projection 
m = N/2+n/A may assume integer, half-integer and quarter-integer values. Thus the wave functions 
of the n-dimensional harmonic oscillator form a basis of the two-valued irreducible representation 
V\_ for the Lie algebra of SU (1, 1) . 

Let us discuss in particular the group-theoretical properties of the wave functions fl2.10p of the 
one-dimensional harmonic oscillator. In view of the definition (11. 3p of the creation and annihilation 
operators, 

ia'^N = TiV^TV-i, -m^^TV = VA^ + l^jv+i, (2.21) 
where we have used the familiar differentiation formulas 

H', (0 = 2kHk^^ (0 = 2^Hk (0 - i^fc+i (0 (2.22) 

for the Hermite polynomials. In this case the relations (12.191) hold for the basis functions of the 
form 

f (-l)^/'^jv, iV = iV+ = 0,2,4, ... forj = -3/4, 

'Ipjm = { (2.23) 
[ N = N- = 1,3,5, ... forj = -l/4, 

where m = N/2 + 1/4. Thus the even and odd wave functions \E'Ar (x) of the one- dimensional 
harmonic oscillator form, respectively, bases for the two irreducible representations P;^ of the algebra 
SU (1, 1) with the moments j = —3/4 for the even values of = A^"*" and j = —1/4 for odd A^ = A^~; 
see ^46j and [55] for more details. 



3. ElGENFUNCTION EXPANSION FOR THE TIME-DEPENDENT SCHRODINGER EQUATION 



In spirit of Dirac's time-dependent perturbation theory in quantum mechanics, see [20], [26], |30j . 
[36] , [S], [12], [53], we are looking for a solution of the initial values problem (11.7l) -( fL9i) as an 
infinite multiple series 

oo 

lP =lP{x,t) = Cm (t) ijjm{u} (X) , (3.1) 

j{u} m=j+l 

where ipjm{u} i^) are the oscillator wave functions (I2.17P depending on the space coordinates x only 
and Cm (t) = Cjm{u} (t) are yet unknown time-dependent coefficients. 

The Hamiltonian (11.81) belongs to a more general type 

H (t) =ujJo + 6 (t) J+ + 6* (t) J_ (3.2) 

with 6 (t) = e"^"^* and uj = 2; see (12.13^ for the definition of operators J± and Jq- It is convenient to 
proceed further with an arbitrary value of the parameter u and then to choose a particular value. 

Substituting expension (13.11) into the Schrodinger wave equation (11.71) . with the help of (I3.2p . 
(I2.19P and the orthogonality property of the oscillator wave functions (I2.17p . namely, 

'^jm{iy} i'j'm'iu'} {x) dv = 6jj'6mm' {S{u}{u'}) , (3.3) 

we obtain an infinite system of the first order ordinary differential equations 

'^^^ = Cm (t) + 6 (t) v/(m-j-l)(m + j) Cm-i (t) (3.4) 
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+ S* {t) v/(m + j + l){m- j) Cm+i (t) . 



The following Ansatz 



Cm (t) = W ] e-'^'^'u^ (t) (3.5) 

results in 

= S (t) e*-* (m + j) M^.i + 6* {t) e~'^' (m - j) (3.6) 
When S (t) = e"*"^* we obtain the system of the first order linear equations 

i— = Au, w(0) = w°, (3.7) 

or, explicitly, 

du 

i-^ = {m + j) Um-i + {m - j) Um+i {m = j + 1, j + 2, ... , oo) (3.8) 

with the Jacobi, or three diagonal, infinite matrix A independent of time t. The system in hand 
should be solved subject to the initial conditions 



< = < (0) = W , c„ (0) (3.9) 

[m — J — 1)1 



in view of ([L9D, and (l33|) . 

The solution of the initial value problem fl3.8l) - fl3.9p can be constructed as follows 

Urn (t) = X^Mm' Umm' (t) , (3.10) 
m' 

where Umm' [t) is a "Green" function, or particular solutions that satisfy the simplest initial condi- 
tions 

Umm' (0) = 5mm'- (3-11) 

Thus the solution of the original initial value problem fll.7p - fll.9l) is given by 

oo 

ip{x,t)=^ ^ Cm{t) lljjm{u}{x) (3.12) 
m=j+l 

with 



(m + jj! 
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where the Green function Umm' (t) will be constructed in this paper in terms of the so-called 
Bargmann function [13], [IS]; see also fl6.8p . It will be done in section 6 after discussion of some 
properties of the Meixner-PoUaczek polynomials and an integral evaluation in the next two sections. 

4. The Meixner and Pollaczek Polynomials 

The exact solution of the initial value problem (13.71) can be obtained with the help of the so-called 
Meixner-Pollaczek polynomials. They can be introduced in the following way. 

The Meixner polynomials [Sn]^[in], [I], [H], [21], [IH] are given by 

where (7),^ = 7 (7 + 1) ■ ... ■ (7 + n — 1) = F (7 + n) /F (7) ; see also [12] for the definition of the 
generalized hypergeometric series. Their orthogonality property is 

f; (k) (k) = -Jf^ S.a (4.2) 

fc=0 ■ H' \ H'J 

with 7 > and < /i < 1; the proof is given, for example, in [16] and |47j . 
An important generating relation for the hypergeometric function 



k=0 



V 



= (1 + sY-''-' (1 + s - suf (1 + s - svY 

^ f -p, -q SUV , ,^ ^, 

X 2Fi\ ; - — — , (4.3) 

\ — r (1 + s — su) (1 + s — sv) 

which is due to Meixner [ID], gives an extension of the orthogonality property as 

f^ml-^'Hk) m^^'Hk) ^-^^^^ (4.4) 

fc=0 

equation (14. 2 p arises in the limit t — 1~, and the explicit representation for the Poisson kernel 

(fit) 

n=0 



^^"^2/) 7^ (4-5) 
w J„ 



(l-/it)"^^^^ y 7 /i(l-t)^ 
The last two equations are related to each other in view of the self-duality 

^i"' ik) I (7). = (n) / (7)fc (n. A: = 0, 1, 2, ... ) (4.6) 

of the Meixner polynomials; cf. (14. ip . 
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The Meixner-PoUaczek polynomials [39], [H], [19], [8], [16] given by 

-in<fi 

p„ (x) = (x, y.) = ^ m(2^' ^) (^x - A) , /i = e~''^ (4.7) 

n! 

with A > and < (p < tt satisfy the following three term recurrence relation 
and the continuous orthogonality relation 

/ {X, <f) P^ (X, ^) p (X) dx = (4.9) 

with respect to the weight function 

p{x) = ^ (2 sin ipf^ |r (A + ix) I' e(2^-")". (4.10) 

The Poisson kernel 

^ wrf = " (r^^ j ^'-''^ 



n=0 

,2 



( X-tx, X-ty _ 4tsii^\ 
X 2F,^ 2A ' -(T^TfJ' 

follows directly from (14.51) and (14. 7p . See [19] for a more general nonsymmetric form of this Poisson 
kernel; its g-expensions are given in P] and [51]. An extension of the orthogonality property is the 
following Fourier integral 

- / e-2^-*P„^(x,v.) P^(x,v.) e(^^-)nr(A + ^x)|^dx (4.12) 
r (2A + n)r (2A + m) e^"^ (sinh t)"+™ 



4^r (2A) n\m\ (cos sinh t + i sin 99 cosh t) 



X 2F1 



-ra, — m / smy^ 



2A ' I 



sinht 



This integral evaluation will be given in the next section. In the limit t — we obtain the 
orthogonality relation (l4.9l) -( HrT0l) . See also [3], [5], [8], and [3l] for the introduction and properties 
of the continuous Hahn polynomials that generalize the Meixner-PoUaczek polynomials. 

5. Evaluation of the Integral 

We consider the analytic continuation of the relation (14.41) in the parameter p. By the Cauchy 
residue theorem the left-hand side of (14. 4p can be rewritten as an integral over the contour Ci (see 
Figure 1), namely, 

-i, / m(7' ^) (z) m(7' {z) e p{z) dz (5.1) 

(i-^r^ p ( -X. -y {^-i^ft 
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Figure 1. Contours in the complex plane. 



where p {z) = (7)^ F (— z) {—pY and |t| < 1. On the semicircle z = — 7/2 + R e*^ with —n/2 < 9 < 
71/ 2 the following estimate 

p{z) = O (i?^-^exp(i?(cos^ln|/i| - sin ^ (arg (-/i) ± tt)))) (5.2) 

holds as -R — s> cx) [8]. Therefore for < 1 and |arg(— yu)| < tt the contour Ci in (15.11) can be 
replaced by the contour C2 where z = —7/2 + ix and — 00 < x < 00. In view of the estimate (15. 2p . 
when |arg (— /i)| < vr the integral in (15.11) converges uniformly on the contour C2, where 6 = ±7r/2, 
for all values of . As a result, this integral can be analytically continued in the parameter /i 
to the entire complex /i-plane with the cut along the positive real axis Re/i > 0. In particular, 
equation (15.11) remains valid for both p = exp {—2iip) and z = —7/2 + ix (7 > 0, < < vr) . The 
substitution (14.71) results in the integral (14.121) evaluation when t — > e^^*. 

6. Solution of the Initial Value Problem 

We can now construct an explicit solution to the original Cauchy problem for the time-dependent 
Schrodinger equation of a modified oscillator in (ll.7p - (ll.9l) . More precisely, we will solve the fol- 
lowing partial differential equation 

^^ = ^E(-(l + ^°^2t) ^ + (l-cos2t) xlA (6.1) 

s=l V / 

subject to the initial condition 



(6.2) 
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by using the eigenfunction expansion fl3.12p - fl3.13p . Looking for a particular solution of the system 
(13. 8p in the form 

u.m (t) = e-'*«* pm (0 , (6.3) 
where ^ is a spectral parameter, one gets 

2ePm (0 = {m+ 3) Pm-1 (0 + {m- 3) Prn+1 (0 (6.4) 

which coincides with the recurrence relation for the Meixner-PoUaczek polynomials (14.81) when 
A = j + 1 and if = tt/2. Thus 

Pm (0 = P^-l f ) • (6.5) 

In view of the orthogonality relation (14. 9p . the Green function Umm' (t) , or solution of the linear 
system (13. 8p that satisfies the initial condition Umm' (0) = 5mm' , can be obtain as the Fourier integral 
over the spectral parameter 

1 f°° 

Umm' it) = :^ e-2*«* Pm (0 Pm' (0 P (0 (6.6) 

"m' "'-oo 

_ (m'-j-l)! 22^+2 
~ {m'+j)\ 27r 

X 



_ (m + j)! (sinht)'"+"'~'^'~' 
- r (2j + 2) (cosht)"^+'"' 

/ -m + j + 1, -m' + j + 1 -L \ 

^ '^'[ 23 + 2 ' -^li^E^J' 

where the last integral has been evaluated with the help of the integral representation f l4.12p . 

The generalized spherical harmonics for the the discrete positive series of the non-compact 
Lorentz group SU (1, 1) are [I3], [16] 

Tim' («, 7) = e— (r) e—'\ (6.7) 
where the Bargmann functions f^^, (r) are given by 

POO 

<m' (^) = {3m\e-'^'^\3m') = e""^/^ / Rnk (r) i?Ari,' (e-^/V) r^-'dr (6.8) 

(_l)'^-i-l / (^^7+7)7(777+7)! /. , r\-2i-2/ ^.m+m' 



r (2j + 2) Y (m - J - 1)! (m' - J - 1)! 
F / -1^ + 3 + 1, -m' +3 + 1 



(smh-j (tanh-j 



2j + 2 ' sinh2(r/2). 

This implies the symmetry relation 

<m'ir) = i-ir-'^'vl,^ir) (6.9) 

and the differentiation formula 

d 



2^^^'m (r) = v/(m-j-l) (m + j) t;^, (r) (6.10) 
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which follows directly from f l6.8p and f l2.19p . Also 

oo 

Yl ^^n,m" (^) ^m'm" (^) = ^mm' , (6.11) 
m"=j+l 

see |13|, HSI, and |56] for more details. 



Combining equations fl3.12p - fl3.13p . 06.60 and 06.80 - 06.90 together, we finally arrive at the follow- 
ing expression in terms of the Bargmann functions 

oo „ 

Cm it) = e-^-™* Yl '"^'"^^'m m / ^ r,m'{u} (^') ^0 ix') dv' (6.12) 

m'=j+l 

for the time-depending coefficients Cm (t) in the expansion fl3.12p . namely, 

oo 

'4){x,t) = Y Y i^jm{u}ix), (6.13) 

m=j+l 

for the solution of the original initial value problem fll.70 - fll.9p with uo = 2. With the help of the 
differentiation formula fl6.10p one can easily verify that these coefficients Cm (t) satisfy the system of 
ordinary differential equations 03. 4p with corresponding initial conditions. This gives a direct proof 
that our solution {x,t) does satisfy the initial value problem fll.7p - fll.90 . 

Thus constructed solution ip {x, t) belongs to the space of the square integrable functions {R"") 
for all times provided that the initial function ipo {x) is of the same class ipo G [R"') ■ Indeed, in 
this case the condition 

oo 

E \^rn{t)f = l (6.14) 
m=j+l 

holds for all values of t in view of the unitary relation 06. lip of Bargmann's functions. We shall take 
advantage of the group-theoretical meaning of this solution in order to cunstruct the corresponding 
Green function later. 



7. The Finite "Rotation" Operators 

We apply a general approach to a certain type of the integral transforms [ST], [SH], [5B], which 
was originated by Wiener [62j, to our problem. Consider the following bilinear sum 

St (r, r') = E (^) (^') (7-1) 

N=K,K+2,... ,oo 

= 2exp (- (r^ + /2) (rr')^ 



oo , . 

2^T(K + n/2 + k) ^"^^^^ 



T{K + n/2 + A;) 
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for the oscillator radial functions {\2.7\i . With the help of the Poisson kernel for the Laguerre 
polynomials 

E (^Ti); ^""^ ^^-^^ 



{x + y)t\ [ - xyt 



see, for example, [50j, p. 212; we derive the following closed form 



bt {r,r ) = -— [rr ) [1 ~ t) exp 



T{K + n/2)' ' ^ 2 l-t 



for this kernel. In the case t = e'" one gets 

C / '\ z, e , 2j+2-n/2 / ' ^ ' 



r(2j + 2)(sin(a/2))'^+' ' ' tan (a/2) 



/ /n2 



where we use the SU (1, 1) moment j = K/2 + n/4 — 1. 
Using the orthogonality property of the radial functions 

RNK{r)RN'K{r) r'^-^dr = 5nn', (7.5) 
from (17. ip we obtain the following integral equation 

tiN-K)/2j^^^ (r) = / St (r, r') Rmk (r') (r')""' rfr', \t\ < 1. (7.6) 
In the SU (1, 1) notations for the oscillator wave functions (12. 171) . when t = e*", it takes the form 

POO 

e'^^'i^.miu} {x) = / Gi (r, r') {x') (r')""' dr' , (7.7) 

where 

2-2i-le*'^(j+l) 2+2-n/2 / ^2_^^/2 



r(2j + 2)(sin(a/2))'^+' ' ' tan (a/2) 



/ /n2 



and x/r = x' /r' = n = n {Q) with = 1. The following symmetry properties hold 

(r, r') = Gi (r', r) = (r, r')) * (7.9) 
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and the formal orthogonality relation is 

{Gi (r, r")yGi (r", r') {r'T' dr" = (7.10) 



r"=0 



where 6 (r) is the Dirac delta function. Also, 

oo 

^'"^"^jmM (x) r.miu} (^') = Yk. {^) Yk. (r, r'), 0. (7.11) 

m=j+l 

In particular, when a = ±vr/2, one gets 

POO 

{±tr ^,m{.} (x) = / Gi,,{r,r')^,^{^y{x') {rT-'dr' (7.12) 

with 

Gi (r,r') = " . , ,±(.^+.-)/2. „fJ , - , ; -1 {rrf] (7.13) 

±7r/2V> ; 2Jr(2j + 2) ^ \^ 2j + 2 ' 2 ^ ^ y ^ ^ 

These formulas will allow us to find a different form of the solution fl6.12p - fl6.13l) in the next section. 

In the process we found out that the finite "rotation" operator e"*""^^ of the group SU (1,1) acts 
on the oscillator wave functions fl2.17l) as the following integral operator 

POD 

e-"''^^ljjm{u} (x) = / Gi, (r, r') (x') (r')""' dr' (7.14) 

Jo 

POO 

= Yk. (n) / Gi, (r, r') R^k {r') {r'T'' dr' 
Jo 

with the kernel explicitly given by (17. 8p . Also, in view of (16. 8p . 

oo 

e-'^'^^jmi,} {x) = J2 (^) ^i^'i-} (^) C^-IS) 

m'=j+l 

see Il6] for more details. 



8. An Integral Form of the Solution 
We rewrite the coefficients f l6.12p in the form 



oo „ 

C„. (t) = e— (-*+'^/^) <'n.m l{{-^r'i^,m'i.}{x')yiJo{x')dv' (8.1) 

m'=j+l -^^^ 

and use the integral transform fl7.12l) as 

poo 



Hr i^.m'iu} ix') = / Gi,, (r', r") ix") {r'T' dr" , (8.2) 

^0 

where x' /r' = x" jr" = n' = n (f2') , in order to obtain 

[t) = e-^™(-*+-/2) (8.3) 
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// poo \ * 

with the help of (17.151) . Substitution into the eigenfunction expansion (16.131) gives 



ij{x,t)=e--'/' iJoix') (8.4) 

yJiui \m=j+l J J 

where by (12X7]) and fTml) 

oo 

E e-^'-^-'+^/'V.M^} (^) r.miu} {e-'x") = Yku m Yk. m GI^,_^/, (r, e-V") , (8.5) 

and our solution takes the form 

^Pix,t)=e-^'/^ [ ^o{x') (8.6) 

JR" 

with uj = 2. The last integral can be evaluated with the help of the formula (115. 5p from the appendix 
at the end of the paper 



Gi^,{r',r") [Gi,^^^,{r,e-V')) (r'T-'dr" (8.7) 

gi7r(j+l) ,^n2j+2-n/2 



^nt/2 



X exp 
X exp 



2^T (2j + 1) (e-* cos (t + 7r/4) - e* sin (t + 7r/4))'^'+' 
e* cos (t + 7r/4) + e~* sin (t + 7r/4) " 



2z e-* cos {t + 7r/4) - e* sin {t + n/A) 

{r'f e~* cos {t + 7r/4) + e* sin (t + 7r/4) 
2i e~*cos(t + 7r/4) -e*sin(t + 7r/4) 

— (rr')^ 



° M 2j + 2 ' 2 (e-* cos (t + 7r/4) - e* sin (t + 7r/4)) 



2 



As a result, the solution of the Cauchy initial value problem (ll.7p - (ll.9l) is given by 

i!{x,t)= / Gt{x,x') ipo{x') dv', 

JR" 



18) 
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where the Green function is 

{x, x') = J2 Yku m Yk. m (r, r') (8.9) 



Ku 

with 



-iTT{2K+n)/A (rr'^^ 



^^'^''''"'^ 2^+W2-ir (ir + n/2) (cos t sinh t + sin t cosh -'^-^""^ ^^'^^^ 
(r^ + (r')^) cost cosh t — (r^ — (r')^) sintsinhi 



X exp i- 



X oFi 



2 (cos t sinh t + sin t cosh t) 

.A 2 



rr 



K + n/2 ' 4 (cost sinh t + sin t cosh t)^ 
The details of the calculations are left to the reader. 

The Green function can also be independently found by separation of the variables in the Cartesian 
coordinates. Indeed, when n = 1 and i^' = 0, 1 with the help of the familiar relations 

„ / — 



cos a = oFi ' — ' sina = a o^i ' — 

our equations flS.Qp -f lHTTOl) can be reduced to 

Gt (x, x') = i {g^ (x, x') + Q\ (x, x')) (8.12) 



2tx% (cos t sinh t + sin t cosh t) 

[x^ — (x')^) sin t sinh t + 2xx' — {xP' + (x')^) cost cosh t 



X exp 



H (cos t sinh t + sin t cosh t) 



which gives the Green function for the one- dimensional Schrodinger equation fl6.ip : we shall elabo- 
rate on the one-dimensional case in the next section. Thus, in the general case, 

n 

) (8.13) 

n/2 



iTii (cos t sinh t + sin t cosh t) 

a;^ — x''^\ sin t sinh t ^2x- x' — (x^ + x'"^) cost cosh t 



X exp , 

2i (cos t sinh t + sin t cosh t) 

and equation ( 18. 9p gives an expansion formula for this Green function in terms of the corresponding 
hyperspherical harmonics. This type of oscillatory integrals is discussed in |57j. 



The time evolution operator for the time-dependent Schrodinger equation (11.71) can formally be 
written as 

U (t, to) = T f exp {-'- I H (f) dA \ , (8.14) 



n 



to 
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where T is the time ordering operator which orders operators with larger times to the left [18], [26]. 
Namely, this unitary operator takes a state at time to to a state at time t, so that 

^{x,t) = U{t,to)^{x,to) (8.15) 

and 

Uit,to) = Uit,t')Uit',to), (8.16) 

U-'{t,to) = UHt,to) = U{to,t). (8.17) 

We have constructed this time evolution operator explicitly in (18.81) . as the integral operator with 
the kernel given by fl8.13p . for the particular form of the time-dependent Hamiltonian operator of a 
modified oscillator in (11. 8p . 

9. The Forced Modified Oscillator 

In the previous section, among other things, we have solved the following one-dimensional time- 
depending Schrodinger equation for a modified oscillator 

on the infinite interval subject to the initial condition 

{x,t)\^^Q = ipo (x) (— oo < a; < oo) . (9.2) 

Our solution takes the form 

/oo 
Go{x,y,t) iJo{y) dy, (9.3) 
-oo 

where the Green function (or Feynman's propagator) is given by 

Go{x,y,t)= ^ = (9.4) 

^y 2ni (cos t sinh t + sin t cosh t) 

(x'^ — y^) sin t sinh t + 2xy — (x^ -|- y^) cos t cosh t 

X exp ' 



2i (cos t sinh t + sin t cosh t) 

This expression may be considered as a generalization of the propagator for the simple harmonic 
oscillator; see [22], 1231, [211, [301, [SD, [HI, and references therein. 



In this section, we shall extend this solution to a more general case of the forced modified oscillator 
with the Schrodinger equation of the form 

^^ = -2(l + cos2t) ^ + -(l-cos2t) x^i, (9.5) 
-^sin2t (^2x|^+ ^-j{t) x^\%g{t) 

where / [t) and g {t) are two arbitrary real valued functions of time only. Indeed, by a method 
similar to one in [2H] and [H] for the case of the forced harmonic oscillator, one can look for the 
Green function in the form 

^ = M e'^, (9.6) 
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where u = Gq (x, y, t) is the fundamental solution of the Schrodinger equation for the modified 
oscillator fl9.1l) and S = a (t) x + (3 (t) y + j (t) . Its substitution into (19.51) results in 

da dB d'-y\ , ^ 

— U 9.7 
dt dt^ dt J ^ ^ 

= - (1 + cos2t) I —a'^u + 2ia—- ] — sin2t axu + fxu + g { au — i—- ] , 
2 \ ox J \ ox J 

where by flOl) 

du .X (costcosht — sintsinht) — 7/ _ ^, 

^=^— : -u- (9.8 

ox cos t smh t + sm t cosh t 

Thus 

da dp d'-y 

— -X H — — w H — — = — sm2tax+fx (9.9) 
dt dt dt ■' ^ ' 

1 , X / 9 X (cos t cosh t — sin t sinh t) — y 
- - 1 + cos 2t) a^ + 2a— — - 

2 \ cos t smh t + sm t cosh t 

/ X (cos t cosh t — sin t sinh t) — y 

+ ^ « + ^ r-; : H ^ 

\ COS t smh t + sm t cosh t 

and equating the coefficients of x, y and 1, we obtain the following system of the first order differ- 
ential equations 

da(t) 2 , , ^ , , / X 1 — tant tanht 

+ . . , . « W = / W + 9 (t) , , , , , (9.10) 



and 



In view of 



dt tan t + tanh t tan t + tanh t 

d/3 (t) _ (1 + cos 2t) a{t)-g (t) 
dt cos t sinh t + sin t cosh t 



(9.11) 



/i'(t) 2 



a(t)c/(t) - -(l + cos2t)a^(t). (9.12) 
2 

/i (t) = cost sinh t + sin t cosh t, (9.13) 



yLi (i) tant + tanht 
equation (19.101) takes the form 

d , , , /x/./s /n1 — tant tanht\ 

- W « W) = , W (/ w + « (i) ) . (9.14) 

The solutions are 

a (if:) = (cos t sinh t + sin t cosh t) (9.15) 

^ / (/ (s) (cos s sinh s + sin s cosh s) + (7 (s) (cos s cosh s — sin s sinh s)) ds, 
Jo 

(3 (t) = f (l + cos2g)a(3) -^(g) 
Jo cos s sinh s + sin s cosh s 

7 (t) = / (s) g{s)-]^{l + cos 2s) (s)") ds. (9.17) 
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The solution of the Schrodinger equation fl9.5p with the initial condition fl9.2p has the form 

/>oo 

^{x,t) = / G{x,y,t)^o{y) dy (9.18) 



with the Green function given by 

G (x, y, t) = Go (x, y, t) e^("(*)-+'3W2/+7W). (9. 19) 

Extensions for the n-dimensional case and to the corresponding heat equation are obvious. The 
details are left to the reader. 

10. Expansion Formula for a Plane Wave 

Equations f l8.9l) -( l8rT0l) and (18.131) imply the familiar expansion formula of a plane wave in HT in 
terms of the hyperspherical harmonics 

n/2 



' (^) " E m Yku m JK+n/2-1 (rrO , (10.1) 



Ku 

where 

is the Bessel function. Although expansion (110. ip is well known in the three dimensional case |26j . 
we were not able to find it in the literature for the (hyper) spherical system of coordinates in i?", 
corresponding to a general binary tree in Vilenkin-Kuznetsov-Smorodinskii's graphical approach 
[16]. Dick Askey has informed us that this formula follows immediately from the Funk-Hecke 
theorem [T] and a familiar integral giving the Bessel functions as a Fourier transform. Bochner 
essentially has this in his book [17]; also see Glaus Miiller's lectures on spherical harmonics [13] 
and [H]; it is probably in some notes of Galderon published in Argentina, but they are not widely 
available. Also see [10], [11], and recent papers [15] and [16] . 

Let us consider a few examples. 

Case n = 1. This is simply Euler's formula 

e'"^ = cos(f + ismip (10.3) 

in view of the familiar relations 

J-i/2{z) = \ — cosz, Ji/2{z) = \ — sinz. (10.4) 

\ TIZ \ TTZ 

Case n = 2. Let x = rcos6 and y = rsin^; see Figure 2(a). Then the expansion formula (110. 11) 
simplifies to 

00 

irr' cosi9-9') _ \^ j (n imid-d') 



or 

giz sin 



^"'Jmirr') e'^^'-''\ (10.5) 
=—00 

00 

Jm {z) e*™'^. (10.6) 
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X y 



m (a) 




Figure 2. Spherical harmonics in i?^ and H? . 

This is a well known relation in the theory of Bessel functions; see, for example, and [17] . 

Case n = 3. If a; = r sin 6' cos*/?, y = r sin 6' sin (yj, z = rcos6 (Figure 2(b)), the expansion formula 
( llO.ip takes the familiar form 



see 



/o \3/2 oo / 

V E E "^'+1/2 (^^0 YL (0, v>) Yim {e\ 



1=0 m=-l 

, and P?] for more details. 



(10.7) 




/' (a) 



Figure 3. Spherical harmonics in R^. 

Case n = 4. Two different system of the spherical coordinates in are 

Xi = r sin if) sin 9 sin 
X2 = r sin %p sin 9 cos ip, 
X3 = r sin ip cos 6*, 
Xi = r cos ij) 

and 

Xi = r sin 9 sin <^i, 
X2 = r sin 6' cos y^i, 
X3 = r cos 6* sin (^2, 
0:4 = r cos6' COSV92, 
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see Figures 3(a)-(b), respectively. The corresponding spherical harmonics are given by 

Ym^ (V-, e,^)=A sin' ^ P//„Y^'' '""'^'^ (cos i^) 

X sinl™l e P^^'jii (cos^) e*™^ (/'>/> , 



(10.8) 



where 



and 



A 



Yi {9,ipi,(p2) 



^{21 + 1) {21' + 2) {I- m)\ {I + m)\ (/' - /)! (/' + / + 1)! 

2'+l™l-2/!r (/' + 3/2) ' 



27r 



iVsinl'^l^ cosl'^lepf^l;|if ./,(cos2^) 



(10.9) 



with 



N 



= l^l|+l^2|, 


l^il 


+ |/2|+2,|/i| + |/2| 


+ 4,... .) 


'(2Z + 2) [(/- 




-|-K2|)/2]![(/H 


h|/i 


l + l^2|)/2]! 


[(^ + l^i| 




l^2|)/2]![(/-|/i| 


+ |/2|)/2]! ' 



respectively. Here Pn"' (^) are the Jacobi polynomials. 
The expansion formulas take the forms 



2 oo I' I 



Z'=0 Z=0 m=-l 



and 



0' (27r)^ 



E E E ^i\i2 Yi hh {0', V'l, ^2) Ji+i irr') , (10.11) 



(;-Ki|-l«2|)/2>0 Z/2>(|Zi| + |«2l)/2>0 



respectively. 




Figure 4. Spherical harmonics in i?". 

The n-dimensional case. The canonical system of hyperspherical coordinates in the Euclidean 
space i?" can be set up as follows 

xi = r sin 9n-i sin 9n-2 ■ ■ ■ sin 62 sin 9i, 
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X2 = r sin 9n-i sin 6'„_2 ■ ■ ■ sin 62 cos 9i, 
X3 = r sin 9n-i sin 6'„_2 ■ • • cos 6^2, 



a;„_i = r sin 6'„_i cos 9n-2, 

with < < 27r, < < vr and k = 2, ... ,n — 1 (Figure 4); see [I] and [IB] for more details. 
The hyperspherical harmonics have the form 

n-l 

Yi^_,i„_,... e^_2, ...,e,) = Al[ sin'-^ P/^^-};^' ^^-^^^^ (e,) (10.12) 

A;=3 

X sinl'^l 02 PlfX^'^' "^'^^"^ (cos ^2) e"'^'' 

{In^l >ln-2> ... >/2> 

where 2jk + 1 = Ik + {k — 1) /2 and A is a normahzing constant. The expansion formula is 

^) E ^'"-^ ^.„-.+n/2-i (rrO (10.13) 

^ ^C-i«n-2... «i (^n-l5 ... ,0i) Yi^_^i^_2_ {9n-l,9n-2, ■■■ , ^l) • 

We leave the details to the reader. 

11. The Schrodinger Equation on the Group SU (1, 1) 

Equation (11. 7p admits the following algebraic generalization. Let J± and Jq be infinitesimal 
operators of the Lie algebra SU (1, 1) and |j, m) = ipjm {x) be a basis of the irreducible representation 
V\_ belonging to the discrete positive series in certain Hilbert space p!3]. We may call the equation 

i^-^ = H{t)^P (ILl) 

with the Hamiltonian of the form 

H{t)=u it) Jo + 6 (t) J+ + 6* (t) J_, (11.2) 

where u [t) is a real valued function of time and 6 {t) is a complex valued function of time, the 
time-dependent Schrodinger equation on the group SU (1, 1) . 

The eigenfunction expansion 

oo 

i){x,t)= E Cm{t) Ipjmix) (11.3) 
m=j+l 

results in the system 

'^^iP' = ^ ^ ^™ + ^ W V(m-j-l)(m + j) Cm-i {t) (11.4) 
+ 6* (t) ^/{m + j + l){m-j) Cm+i (t) 
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with m = j + 1, j + 2, ... similar to (13.41) . When 6 (t) = e~"^* and uj =constant, the same consider- 
ations as in sections 3 and 6 give us the solution of the corresponding Cauchy initial value problem 
in the form of the expansion (111.31) with the time-dependent coefficients 



oo 



c™(t)=e--* J2 ^"''-"'vL'mm (^,w, ^o), (11.5) 

m'=j+l 

where f (2t) are the Bargmann functions (16. 8p and 



{if,x)= ^*{x)x{x) dfi{x) (11.6) 

J Supp fi 

is the inner product in the Hilbert space under consideration. A formal substitution of (Ill.Sp into 
(I11.3P results in 



i/j{x,t)= G{x,y,t) ipoiy) dn{y) (11.7) 

J Supp fM 

with 

oo oo 

G{x,y,t)= ^ e-'^^U"^'-^vl,^{2t) ^,rn{x)^,„,,{y). (11.8) 

m=j+l m'=j+l 

The original Cauchy initial value problem (ll.7p - (ll.9p provides an explicit model of the abstract 
Hilbert space. Another realization is a quasipotential model of the relativistic oscillator; see, for 
example, [37j, [7], and [8]. We shall discuss this model and its generalization in the next two 
sections. Propagators for the difference models of the simple harmonic oscillator [9] can be derived 
in a similar fashion. The details are left to the reader. 

12. Propagator for the Relativistic Oscillator 

For the consistent three-dimensional description of a relativistic two-particle system in quantum 
field theory a quasipotential approach had been formulated and, in the framework of this approach, 
some relativistic generalizations of the exactly solvable problems of quantum mechanics have been 
considered; see |37j, [7], and references therein for more information. We shall discuss here a model 
of the (simple) relativistic oscillator described by the following Hamiltonian operator 

Ho = nic^ cosh (iXdx) + -mu'^x (x -|- iX) exp (iXdx) , (12-1) 

where A = H/mc is the Compton wave length, = d/dx and exp (ad^) f (x) = f {x + a) is the 
shift operator. The square-integrable solutions of the stationary Schrodinger equation 

Ho^Hn (x) = S„^n (x) (-00 < o: < oo) (12.2) 

on the real line, which correspond to the discrete energy levels 



1 



E„ = Huj(n + i^), y = - + ^-+(^— j (n = 0,l,2,... ), (12.3) 
can be found in terms of the special Meixner-PoUaczek polynomials [8] as follows 



^n(x)=2^/^-^^---— (z.(z.-l))--/^^ r{u + zx/X) P:{x/X, rr/2). (12.4) 
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On the other hand, solution of the Cauchy initial value problem for the time-dependent Schrodinger 
equation 

iH^ = Ho^ (12.5) 
with the Hamiltonian (112.11) subject to the initial condition 

V'(a;,t)|,=o = ^o(a;) G £'(-oo,oo) (12.6) 

has the form 

OO «QQ 

ij {x, t) = Y, e"^^""" (^) / K (y) ^0 (y) dy (12.7) 

n J -co 



n=0 



Go{x,y,t) ipoiy) dy, 
where the Green function (or Feynman's propagator) 

OO 

Go (x, y,t) = J2 e-*-('^+^)* vl/„ (a;) v]/; (y) (12.8) 

n=0 

can be found in a closed form as follows 

Go{x,y,t) = ^ {^Y''^^^' (^sin(.;t/2))^(^--)/^(cos(a;t/2))-^(^+^^/^ (12.9) 

X r(.(a:-y)/A) ,F,[^ lU{y-x)IX ' ^^^Z^) 
+ -\ ( ^Y""'^'^ sin (cut/2) )*("-^)/^ (cos (cut/2) )-'("+^)/^ 
X Tiiiv-xMX) r('^ + WA)r(z/-^|//A) 

^ F ( ^-^yl^^ l-u-iy/X _ -^2/ 

We have used here the Poisson kernel for the Meixner-Pollaczek polynomials (14.111) in order to sum 
the series and a transformation formula for the analytic continuation of the hypergeometric function 
[I] and [Sj. 

13. A Modified Relativistic Oscillator 

As is known fT] a dynamical symmetry group for the relativistic oscillator with the Hamiltonian 
(112. ip is the group SU (1,1) (or isomorphic groups S'O (2, 1) ~ Sp{2,R) ~ SL {2, R)), whose 
generators are realized as the difference operators 

1 T 1 C 

Ko = —Ho, K± = -±iKoT^exp{-iXd^) (13.1) 
fiw X Xuj 

acting on the eigenfunctions of the Hamiltonian (112. ip as follows 



K+^Hn (x) = ^(n + 1) (n + 2z/) (x) , (13.2) 

(x) = y/n {n + 2iy- 1) (x) , 

i^o^n (x) = {n + iy) ^„ (x) . 
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The wavefunctions (x) = ipjm with n = m — j — 1 and u = j + 1 form the basis for the infinite- 
dimensional irreducible unitary representations of the discrete positive series P;^ of the universal 
covering group SU (1,1); see equations f l2.19p . 

In this section, we shall solve the time-dependent Schrodinger equation 

tn^=H{t)^P (13.3) 

with a modified Hamiltonian of the form 

H{t) = Ho + h {6 (t) K+ + 6* (t) K_), 6{t) = e"^"* (13.4) 
subject to the initial condition (112.61) . The eigenfunction expansion 

oo 

^{X,t)= Cm{t) ^jm{x) (13.5) 

m=j+l 

leads to the familiar system (13.41) . whose solutions in terms of the Bargmann functions are 

^ POO 

c„(t)=e— * Yl '"'''"'^'n.m / r,m'iy)i^oiy) dy. (13.6) 

m'=j+l 

Thus 

/oo 
G{x,y,t) iJo{y) dy, (13.7) 
-oo 

where the Green function is given by 

oo oo 

G{x,y,t)= ^"''"''^'"''^"'^'n.m ^,^{x)r,m'iy)- (13.8) 

m=j+l m'=j+l 

This multiple series can be simplified to two single sums. Indeed, rewriting the wave functions 
(112. 4p in terms of the hypergeometric function 



2^'^' (j(j + l))-/^^r(j + l + .x/A) / (m + j)! 
by (14.71) and (14. ip and then using the Meixner generating relation (14. 3p . we have 

oo 
m'=j+l 



V2^ r {2j + 2) 



e 



^^"'^ V(--J-l)! ^ 2, + 2 
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Thus, say for a pure imaginary time, 



oo 



G(a:,y,t)= {-^r~'~' 4^,m{x) (13.11) 



m=j+l 



oo 



m' =j+l 

1 ^. ^ ^-^y{y-x)/2\ r (j + 1 + ix/ A) T (j + 1 - i?//A) 



27rA ^-^ ^-^ r (2j + 2) 

X e-^''yl^ (cos (cut/2))-'^-' (« tan {ut/2) fy-''^'^ 
i + 1 — ix/\, j + 1 + iy/\ _ 1 



X 2F1 



23 + 2 ' cos2 (cut/2) 



Here we use a transformation formula for the analytic continuation of the hypergeometric function 
[T] and [3^ in order to obtain the final result as follows 

G (x, y,t) = -^( ^y^'~"^^^ e-^ity/x ^.^ (u;t/2))^(^-^)/^ (cos (cut/2))-'("+^)/^ (13.12) 
27rA WujJ 

X r (. (X - ,) /A) ^1 _ ; sm^ {ut/2) 

+ ^ f-^V^'"''^^' (2sin(cut/2))^("-^)/" (cos (cut/2) )*("+^)/' 

27rA WuJ/ 

xTH(v-x) /X) r(^ + WA)r(^-zi//A) 

/ u + ix/X, l-u + ix/X . 2/ ,/o\ 
^ ^^^(, l + .(x-y)/A ' 

Comparing expressions fll2.9l) and (113.121) for two propagators, one can see that 

G{x,y,t) = Go{x,y,t) e-2^*^/\ (13.13) 

It also follows from mBf . ^M . ffT3J0D . and ffmi]) . The details are left to the reader. 



14. A System of Ordinary Differential Equations 
Consider the following infinite system of ordinary differential equations 

= (^ + 1) Un+i (t) -2{n + X) cos(^ Un (t) + (n + 2A - 1) u^-i [t) , (14.1) 

with t>_i (t) = (ra = 0, 1, 2, ... ) subject to the initial conditions 

UniO) = ul (14.2) 
Looking for a particular solution in the form 

un{t) = e-'"'^'-^ P::{x,^), (14.3) 
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one gets the three term recurrence relation (14 .Sp for the Meixner-PoUaczek polynomials (x, (p) . 
The particular solution that satisfies the initial conditions Vnm (0) = 5nm is given in the integral 
form 

1 /"°° 

Un^ it) = -^ {x, if) Pi (x, if) p (x) dx, (14.4) 

where p(x)and c?^ are the weight function and the squared norm for the Meixner-Pollaczek poly- 
nomials P^ (x, (f) , respectively; see (I4.9p -( 14T0|) . The solution of the initial value problem is 



(t) = 5^ < U^m {t) . (14.5) 



m=0 

The last integral can be evaluated as a single sum with the help of Meixner's generating relation 
(14. 3p . The final result is 

(2A)„ e-Msiny.)^^ (sinh (t sin y))-+- 

(cos (fi smh [t sm (p) + t sm (p cosh (r sm (p) ) 

^ ^ / —n, —m _ / sin(y9 ^ ^ 
^ ^ I 2 A ' ysinh (t sin (/?) 

This is an extention of (16. 30 - 06. 6p where ip = 7r/2. 
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15. Appendix. Another Integral Evaluation 
The following integral 

iFJ " ; kzj iFi ( ; k'zj dz (15.1) 



r (7) A"+"'-^ (A - A;)"" (A - 2F1 



a, a 



kk' 



7 ' {\-k){\-k') 

is evaluated in [29] and [36]. When a' = a, replace A = fia, z = x/a and take the limit a = n — > cxo 
with the help of 

hm 1^1 f " ; -) = oi^i f : ; kx) , (15.2) 
_ / (T. (T kk' \ _ / — kk'\ , . 



n— +0O \ 77, 

The end result 
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lim ( 1 + = e^. (15.4) 



= r (7) y^-Te(*=+'=')/^ oi^i 



- _ w_ 

7 ' 7^ 



is required in section 8 of this paper. 
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